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Abstract 

The Rutherford Backscattering Spectroscopy (RBS) and Particle Induced X-ray 
Emission (PIXE) techniques can be used to obtain complementary information 
about the characteristics of a sample but, traditionally, a gap has separated the 
available computer codes for analyzing data from each technique, being hard to si- 
multaneously analyze data from the same sample. The recent development of a free 
and open source library, LibCPIXE, for PIXE simulation and analysis of arbitrary 
multilayered samples, has permitted to integrate this technique into the DataFur- 
nace code which already handles many other IBA techniques such as Rutherford 
and non-Rutherford backscattering, elastic recoil detection, and non-resonant nu- 
clear reaction analysis. The fitting capabilities of DataFurnace can therefore now 
be applied to PIXE spectra as well, including the Bayesian Inference analysis and 
the simultaneous and coherent fitting of multiple spectra from different techniques. 
Various examples are presented in which the simultaneous RBS and PIXE analy- 
sis allows us to obtain consistent results that cannot be obtained by independent 
analysis of the data from each technique. 
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1 Introduction 



Despite of the close relation of particle-based IBA techniques (RBS, ERDA,. . . ) 
with photon-based ones (PIXE, PIGE,. . . ), and despite the fact that they are 

often used to obtain complementary information about the same sample, the 
most popular IBA analysis codes are focused either on one or the other cate- 
gory and data cannot be analyzed in a coherent and unified way. An attempt 
to unify PIXE and RBS analysis had been already done [1] by building a 
common interface to the well known RUMP [2] and GUPIX [3] programs. In 
this work we present a more ambitious initiative in terms of integrated anal- 
ysis of IBA techniques: the implementation of the support for FIXE in the 
DataFurnace code (which already supported RBS, ERDA, and non-resonant 
NRA), achieved by using the recently developed LibCPIXE free hbrary for 
PIXE yield simulation [4] . This imphes that all the fitting capabihties already 
available in DataFurnace for other techniques can now be applied to FIXE, 
including simulated annealing fits and the Bayesian inference error estimation. 



2 DataFurnace + LibCPIXE 



The integration of LibCFlXE with DataFurnace, incorporates FIXE to the 
repertory of supported techniques in the latter code. This opens up the possi- 
bility of simultaneously analyzing data from several spectra of a given sample 
and thus merging the information from all of them, producing a self-consistent 
result. Note that the experimental spectra can be obtained with different tech- 
niques (now including FIXE) as well as with the same technique under different 
experimental conditions in order to obtain complementary information about 
the sample. 

In contrast with the most widely used programs for FIXE data analysis, which 
make use of nonlinear least squares or Marquardt fitting algorithms [5], the 
DataFurnace code implements simulated annealing and Bayesian inference al- 
gorithms. The combination of these algorithms has proven to be a very pow- 
erful tool when dealing with IBA problems in which one may have multiple 
and/or suboptimal formal solutions [6]. In relation with FIXE analysis, this 
provides a very flexible fitting procedure that, e.g, does not pose any limit to 
the analysis of layered samples with a given element present in more than one 
layer (which is, currently, a limitation in the GUPIX code [7]). The Bayesian 
inference has also the advantage of providing a mathematically rigorous esti- 
mation of the errors in the fitted parameters (concentration profiles for each 
element in the sample). 

The treatment of the PIXE data in DataFurnace differs from that done for 
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other techniques or that done by other PIXE analysis codes: instead of fitting 
the PIXE spectra channel-by-channel, only user-selected characteristic X-ray 
peak areas are fitted (avoiding the difficulties of fitting background and/or 
detector response functions and allowing to discard peaks which are non-apt 
for the fit). In order to identify the peaks and filter the background, the 
QXAS program [8] can be used (DataFurnace can directly read its output files) 
although a more convenient program is currently being developed and will be 
incorporated into the LibCPIXE library. Besides this detail, the treatment 
of the PIXE spectra in DataFurnace is identical to that done for the other 
techniques. See [9] for a detailed description on how the fitting process works 
in DataFurnace, and how the information provided by different techniques 
about a given sample is merged. 

A typical case of ambiguity found in PIXE spectra is related to the occur- 
rence of non-resolved overlapping peaks. This can be handled either by simply 
ignoring them in the fit (and relaying on other peaks for those elements) or 
by fitting the sum instead of the separated values. Note that in many cases, 
the ambiguities can be also ruled out by providing additional spectra (PIXE 
and/or other techniques) or by setting constrains to the fit based on known 
information about the sample. 



3 Examples 

A first example consists in obtaining the concentration of the elements present 
in a bulk stoichiometric Sb2S3 sample. A PIXE spectrum with 1.1 MeV H+ 
was acquired using = 15° and (pdet = 55°, where and ^^^et are, re- 
spectively, the incidence and detection angles relative to the surface normal. 
A Bayesian fit of the areas of the Sb i^a^ 2 and S Ka^ .^ peaks was performed 
(the only free parameters in this case being the bulk concentration of S and 
Sb as well as the beam flucnce). The results, as shown in Table 1, reproduce 
the expected stoichiometric values and provide a reasonable error estimation 
refiecting mainly the statistical uncertainties in the experimental data. 

A second example consists in the analysis of a layered Mn4lr/Si sample (a 
Mn4lr film on top of a Si substrate). A Bayesian fit of the Mn Ka^ 2) 1^ -^ai 2 
and Si lines of a 1.22MeV PIXE spectrum (with (pine = 7.5° and 

(pdet = 77.5°) was performed leaving the thickness and composition of the film 
as free parameters. The results are shown in Table 1 and are compatible with 
the nominal values (stoichiometric concentrations and a previous RBS thick- 
ness measurement) within the errors estimated by the Bayesian inference. In 
order to ascertain the accuracy of these results, we performed a simultaneous 
RBSh-PIXE analysis using the same PIXE data as before and a 2 MeV He"'" 
RBS spectrum taken at normal incidence angle and 177° scattering angle. As 
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expected, the uncertainties are greatly reduced in this case (see Table 1) since 
the RBS data introduces a tight constrain in the thickness of the layer. 

When dealing with a more complex problem such as obtaining the concen- 
tration profiles for all elements in a GalnAsSb/GaSb sample, the informa- 
tion from a single PIXE spectrum (or a single RBS spectrum) is not enough 
to solve the ambiguities of the fit. Previous analysis of this kind of samples 
[10,11] involved the separated analysis of RBS and PIXE spectra and a manual 
iteration to obtain coherent results from both techniques. The DataFurnace 
code now allows to perform the analysis in a much more convenient way: in 
the third example, the simultaneous analysis of six PIXE spectra taken with 
different beam energies and incidence angles proves useful for characterizing 
the sample. The Bayesian fit (shown in Fig. 1) was performed assuming a 
Gaa;Ini_a,Asi_ySbj^ {x,y < 1) of unknown thickness and composition on top of 
a stoichiometric GaSb substrate. The results, shown in Table 1, are compatible 
with those of previous characterizations of the same sample (as seen in Table 

1)- 

Finally, another Bayesian fit has been performed incorporating a 2 McV 
RBS spectrum (see Fig. 2) to the mentioned set of six PIXE spectra. The 
results in Table 1 show that less uncertainties are present in this case compared 
to the PIXE-only one due to the complementarity of the information provided 
by each technique: PIXE is good at determining the amount of each clement 
but has poor depth sensitivity while RBS provides good depth resolution but 
cannot separate elements of similar mass (such as Ga-As or In-Sb). It is 
interesting to remark the fact that this analysis has been performed in one 
single step, providing all the experimental data to the DataFurnace code and 
obtaining a self-consistent result which is more reliable and considerably less 
time consuming than that obtained by using separated analysis codes. 



4 Conclusions 



We have shown that the integration of LibCPIXE into the DataFurnace code 
facilitates very interesting analysis modes such as the simultaneous PIXE -t- 
RBS or the multiple (differential) PIXE analysis. These combined analysis 
modes are, in many cases, the only way of obtaining a consistent characteri- 
zation of a specimen. The use of robust fitting — simulated annealing — and 
rigorous error estimation routines — Bayesian inference — makes it possible to 
characterize arbitrary layered samples, with no limitations on the composition 
of each layer. 
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List of Figures 



1 Comparison between experimental and fitted main lines of the 
four elements in the film. E is the proton beam energy. In 
geometry Gl, (pine = 7.5° and (pdet = 62.5° whereas in geometry 
G2, (pine = 22.5° and (pdet = 47.5°. Spectra obtained using a 
proton beam of 1.05 MeV were detected without any special 
filter between the chamber window and the Be window of the 
detector. The remaining spectra were measured using a 1mm 
thick Mylar filter in front of the x-ray detector. 

2 Experimental RBS spectrum (squares) for 2 MeV H"*" on the 
same GalnAsSb/GaSb sample as in Fig. 1. The fit (sohd line) 
was performed considering not only the RBS spectrum but 
also the data from six FIXE spectra (those of Fig. 1). Partial 
spectra for each element (non-solid lines) are also shown. 
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Fig. 1. Comparison between experimental and fitted main lines of the four ele- 
ments in the film. E is the proton beam energy. In geometry Gl, (pi^c = 7.5° and 
4>det = 62.5° whereas in geometry G2, (pine = 22.5° and ^^et = 47.5°. Spectra ob- 
tained using a proton beam of 1.05 MeV were detected without any special filter 
between the chamber window and the Be window of the detector. The remaining 
spectra were measured using a 1mm thick Mylar filter in front of the x-ray detector. 
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Fig. 2. Experimental RBS spectrum (squares) for 2 MeV on the same 
GalnAsSb/GaSb sample as in Fig. 1. The fit (solid line) was performed consid- 
ering not only the RBS spectrum but also the data from six PIXE spectra (those 
of Fig. 1). Partial spectra for each element (non-solid lines) are also shown. 
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1 Results for the characterization of the examples of section 
3. The "nominal" values are the stoichiometric ones (when 
known) or had been obtained in previous independent 
characterizations of the same samples (marked with '*'). 
The '(S)' superindex indicates "Substrate" (with known 
stoichiometry) . |lO 
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Tabic 1 

Results for the characterization of the examples of section 3. The "nominal" values 
are the stoichiometric ones (when known) or had been obtained in previous indepen- 
dent characterizations of the same samples (marked with '*'). The '(S)' super index 
indicates "Substrate" (with known stoichiometry). 
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